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The spectral theory of tensors is an important part of numerical multi-linear algebra, 
or tensor computation ESI 90J. 

It is possible that the ideas of eigenvalues of tensors had been raised earlier. However, 
it was in 2005, the papers of Lim and Qi initiated the rapid developments of the spectral 
theory of tensors. In 2005, independently, Lim [57] and Qi (70] defined eigenvalues and 
eigenvectors of a real symmetric tensor, and explored their practical application in deter- 
mining positive definiteness of an even degree multivariate form. This work extended the 
classical concept of eigenvalues of square matrices, forms an important part of numerical 
multi-linear algebra, and has found applications or links with automatic control, statisti- 
cal data analysis, optimization, magnetic resonance imaging, solid mechanics, quantum 
physics, higher order Markov chains, spectral hypergraph theory, Finsler geometry, etc, 
and attracted attention of mathematicians from different disciplines. After six years' 
developments, we may classify the spectral theory of tensors to 18 research topics. Be- 
fore describing these 18 research topics in four groups, we state the basic definitions and 
properties of eigenvalues of tensors. 

A. Basic Theory 

An 72,-dimensional homogeneous polynomial form of degree m, f(x), where x G 3ft™, 
is equivalent to the tensor product of a symmetric n-dimensional tensor A = (a^....^) 
of order m, and the rank-one tensor x m : 

n 

f(x) = Ax .= ^ ] a h,— ,im X h ' ' ' x im ■ 

The tensor A is called symmetric as its entries a^...^ are invariant under any per- 
mutation of their indices. The tensor A is called positive definite (semidefinite) if 
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f(x) > (f(x) > 0) for all x G 9ft™, x 7^ 0. When m is even, the positive definiteness of 
such a homogeneous polynomial form f(x) plays an important role in the stability study 
of nonlinear autonomous systems via Liapunov's direct method in Automatic Control 
[66J. For n big and m > 4, this issue is a hard problem in mathematics. 

In 2005, Qi [70] defined eigenvalues and eigenvectors of a real symmetric tensor, 
and explored their practical applications in determining positive definiteness of an even 
degree multivariate form. 

By the tensor product, Ax" 1-1 for a vector x G 9ft n denotes a vector in 3ft n , whose ith 
component is 

n 

(AX )j = ^ ] a i,l2,--- ,im X i2 ' ' ' X im 1 

We call a number A G C an eigenvalue of *4 if it and a nonzero vector x G C n are 
solutions of the following homogeneous polynomial equation: 

[Ax™- 1 ) . = Xxf- 1 , Vi = l,...,n. (1) 

and call the solution x an eigenvector of A associated with the eigenvalue A. We call 
an eigenvalue of A an H-eigenvalue of A if it has a real eigenvector x. An eigenvalue 
which is not an H-eigenvalue is called an N-eigenvalue. A real eigenvector associated 
with an H-eigenvalue is called an H-eigenvector. 

The resultant of is a one-dimensional polynomial of A. We call it the charac- 
teristic polynomial of A. 

We have the following conclusions on eigenvalues of an mth order n-dimensional 
symmetric tensor A: 

Theorem 1. (Qi 2005) 

(a) . A number A G C is an eigenvalue of A if and only if it is a root of the charac- 
teristic polynomial (p. 

(b) . The number of eigenvalues of A is d = n{m — l) n_1 . Their product is equal to 
det(A), the resultant of Ax m ~ l = 0. 

(c) . The sum of all the eigenvalues of A is 

{m-l) n -hx{A), 

where tr(A) denotes the sum of the diagonal elements of A. 

(d) . If m is even, then A always has H- eigenvalues. A is positive definite (positive 
semidefinite) if and only if all of its H-eigenvalues are positive (nonnegative). 

(e) . The eigenvalues of A lie in the following n disks: 

|A - a M) ...,i| < 5^{|a iii2 ,... )im | : i 2 ,- ■ ■ ,i m = 1, • • • ,n, {i 2 , ■■■ ,i m } ^ {i, ■ ■ ■ ,i}}, 
for i = 1, ■ • • , n. 
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The resultant of Ax" 1 ' 1 = 0, denoted as det(*4.) in this theorem, is called the sym- 
metric hyperdeterminant of A. When m > 3, it is different from the hyperdeterminant, 
introduced by Cayley [12] in 1845, and studied by Gelfand, Kapranov and Zelevinsky [29] 
in 1994. When m = 2, both the hyperdeterminant and the symmetric hyperdeterminant 
are the same as the determinant in the classical sense. 

This theorem extended the classic theory of matrix eigenvalue to tensors. However, 
eigenvalues and H-eigenvalues are not invariant under orthogonal transformation. This 
is needed in mechanics and physics. Hence, in [70J, another kind of eigenvalues of tensors 
was defined. 

Suppose that A is an mth order n-dimensional symmetric tensor. We say a complex 
number A is an E-eigenvalue of A if there exists a complex vector x such that 

r Ax" 1 - 1 = \x, 

\ x T x = 1. 

In this case, we say that x is an E-eigenvector of the tensor A associated with the E- 
eigenvalue A. If an E-eigenvalue has a real E-eigenvector, then we call it a Z-eigenvalue 
and call the real E-eigenvector a Z-eigenvector. 
When m is even, the resultant of 

Ax" 1 ' 1 - X{x T x) m ^ 1 x = 

is a one dimensional polynomial of A and is called the E-characteristic polynomial of 

A. We say that A is regular if the following system has no nonzero complex solutions: 

r Ax" 1 - 1 = o, 

\ x T x = 0. 

Let P = (pij) be an n x n real matrix. Define B = P m A as another mth order 
n-dimensional tensor with entries 

n 

^il,»2,— Am ~ ^ ] PhjlPi2j2 ' ' ' Pimjm a jl,j2,--- ,jm- 

jl,j2,— Jm = l 

If P is an orthogonal matrix, then we say that A and B are orthogonally similar. 

Theorem 2. (Qi 2005) 

We have the following conclusions on E-eigenvalues of an mth order n-dimensional 
symmetric tensor A: 

(a) . When A is regular, a complex number is an E-eigenvalue of A if and only if it 
is a root of its E-characteristic polynomial. 

(b) . Z- eigenvalues always exist. An even order symmetric tensor is positive definite 
if and only if all of its Z- eigenvalues are positive. 
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(c) . If A and B are orthogonally similar, then they have the same E-eigenvalues and 
Z-eigenvalues. 

(d) . If A is the Z-eigenvalue of A with the largest absolute value and x is a Z- 
eigenvector associated with it, then Xx m is the best rank-one approximation of A, i.e., 

\\A - Xx m \\ F = ^H^lll^ - A 2 = min{||^ - au rn \\ F : a G % u G 5T, ||u|| 2 = 1}, 

where || • \\p is the Frobenius norm. 

Theorem[2](d) indicates that Z-eigenvalues play an important role in the best rank-one 
approximation. This also implies that Z-eigenvalues are significant in practice. The best 
rank-one approximation of higher order tensors has extensive engineering and statistical 
applications, such as Statistical Data Analysis [2"5 ] W7 \ [TUB] . 

Eigenvalues, H-eigenvalues, E-eigenvalues and Z-eigenvalues were extended to non- 
symmetric tensors in [72]. When m is odd, the E-characteristic polynomial was defined 
in [72] as the resultant of 

( A T m-l _ \ m-2 

j nn~ T rr> /y»2 

Independently, Lim also defined eigenvalues for tensors in [57]. Lim defined eigen- 
values for general real tensors in the real field. The I 2 eigenvalues of tensors defined 
by Lim are Z-eigenvalues of Qi [70], while the l k eigenvalues of tensors defined by Lim 
are H-eigenvalues in Qi [70]. Notably, Lim proposed a multilinear generalization of the 
Perron-Frobenius theorem based upon the notion of l k eigenvalues (H-eigenvalues) of 
tensors. 

The calculation of eigenvalues is NP-hard with respect to n when m > 3. In the case 
of Z-eigenvalues, when n is small, we may solve polynomial system (J2J) by elimination 
methods, such as the Grobner method. When n = 2 or 3, we may eliminate A in 
Ax m ~ l = Xx first. See [77J. When n is not small, we may use the power method. This 
method only can get one eigenvalue and cannot guarantee it is the largest or the smallest 
eigenvalue. See [49J. 

B. Eigenvalues of Nonnegative Tensors, Higher-Order Markov Chains and 
Spectral Hypergraph Theory 

This group include four research topics. 

(1) The Largest H-Eigenvalue of a Nonnegative Tensor. We call p(A) the 
spectral radius of tensor A if 

p(A) = max{|A| : A is an eigenvalue of A}, 

where |A| denotes the modulus of A. A tensor A = (oj r --i m ) is called nonnegative, if 
a h-i m ^ 0, for all i\, . . . ,i m = 1, • • • ,n. In the study of the Perron-Frobenius theorem 
for nonnegative tensors and related algorithms, several classes of (nonnegative) tensors 



4 



were generalized to (nonnegative) tensors, some new classes of (nonnegative) tensors 
were introduced. Among these, Chang, Pearson and Zhang [15] extended irreducible 
matrices to irreducible tensors; Pearson [69J extended essentially positive matrices to 
essentially positive tensors; Chang, Pearson and Zhang [16] extended primitive matrices 
to primitive tensors; Friedland, Gaubert and Han [28] introduced weakly irreducible 
tensors and wekly primitive tensors; Zhang and Qi |102j introduced weakly positive 
tensors; Hu, Huang and Qi [38] introduced strictly nonnegative tensors and described 
the relations among corresponding nonnegative tensors. 

Now we may state the Perron- Frobenius theorem for nonnegative tensors as follows. 

Theorem 3. (The Perron-Frobenius Theorem for Nonnegative Tensors) If A 

is a nonnegative tensor of order m and dimension n, then p(A) is an eigenvalue of A 
with a nonnegative eigenvector x G 9?" . (Yang and Yang [96]) 

If furthermore A is strictly nonnegative, then p{A) > 0. (Hu, Huang and Qi J3%l) 

If furthermore A is weakly irreducible, then p(A) is an eigenvalue of A with a positive 
eigenvector x G (Friedland, Gaubert and Han 12^ ) 

Suppose that furthermore A is irreducible. If A is an eigenvalue with a nonnegative 
eigenvector, then A = p{A). (Chang, Pearson and Zhang fTh^ ) 

In this case, if there are k distinct eigenvalues of modulus p{A), then the eigenvalues 
are p(A)e l27T ^ k , where j — 0, • • • , k — 1. The number k is called the cyclic index of A. 
(Yang and Yang f9b^) 

If moreover A is primitive, then its cyclic number is 1. (Chang, Pearson and Zhang 

m) 

If A is essentially positive or even-order irreducible, then the unique positive eigen- 
value is real geometrically simple. (Pearson 2010) (Yang and Yang \9J$ ) 

Chang, Pearson and Zhang [15] extended the well-known Collatz minimax theorem 
for irreducible nonnegative matrices to irreducible nonnegative tensors. Based upon 
this, Ng, Qi and Zhou proposed an algorithm for calculating p(A) in [61]. For a weakly 
irreducible nonnegative tensor A, this algorithm generates two sequence {A fc } and {A^} 
such that 

Ai < A 2 < • • • < A = p(A) < • • • < A 2 < Ai. 

The convergence of the NQZ algorithm was established by Chang, Pearson and Zhang 
[T6] under primitivity, and by Friedland, Gaubert and Han [28] under weak primitivity. 
The linear convergence of the NQZ algorithm was established by Zhang and Qi [102J. 
Further developments in this topic can be found in [HI ED EU E3 EB3 S3 EUl HDSl 

fTMfmllml fm]. 

(2) Higher-Order Markov Chains and and Z-eigenvalues of Nonnegative 
Tensors. Ng, Qi and Zhou put higher-order Markov chains as an application in their 
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paper. A higher-order Markov chain model is used to fit the observed data through the 
calculation of higher-order transition probabilities: 



< Pkik 2 -k m = Prob(X t = ki | X t -x = k 2l . . . , X t ^ m+1 = k m ) < 1 (3) 

where 

n 

^2Pk u k 2 ,-,km = !■ ( 4 ) 
fei=i 

Thus, we have an m-order n-dimensional nonnegative tensor V consisting of n m 
entries in between and 1: 

V = (Pfcifca-fem)) 1 - fc l> fc 2, k m < U. 

Let the probability distribution at time t be E S n , where S n = {x e 3fJ n : x > 
0, Yjj=x x i = ■*-}• Then we have 

x m = Vx *-d . . . x ( t - m+ i) := f Pik2 ... km xi; l) ■ ■ ■ 4r +1) ) e ^ (5) 

for t = m + 1 , ■ ■ ■ . Assume that 

lim x {t) = x*. (6) 

t— >oo 

Then we have x* G S'n and we may call x as the stationary probability distribution of 
the higher-order Markov chain. We see x* satisfies the following tensor equation: 

Vx m ~ l = xeS n . (7) 

Li and Ng j55] showed that if V is irreducible, then a positive solution of (I7j) exists. 
Then also give sufficient conditions for uniqueness of the solution of the fixed point 
problem (J7J). Under the same conditions, Li and Ng [56J established linear convergence 
of the power method 

y ^)=V(yM) m ~ l (8) 
Li and Ng [56J established linear convergence of the power method (jSJ), but has not 
proved the convergence of the higher order Markov chain ([5]), i.e., fl6]). For second-order 
Markov chain, under some conditions, Hu and Qi |JT] established the convergence ()6]). 

The fixed point problem (J2J) is different from the eigenvalue problem ([T|). It is some- 
what connected with the Z-eigenvalue problem of a special class of nonnegative tensors. 
Unlike the largest H-eigenvalue of a nonnegative tensor, the largest Z-eigenvalues of a 
nonnegative tensor are not unique in general [15J. The recent paper of Chang, Pearson 
and Zhang [17] reveals some similarities as well as differences between the Z-eigenvalues 
and H-eigenvalues of a nonnegative tensor. 
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(3) Spectral Hypergraph Theory. Lim [57] pointed out that a potential appli- 
cation of the largest eigenvalue of a nonnegative tensor is on hypergraphs. A graph 
can be described by its adjacency matrix. The properties of eigenvalues of the adjacency 
matrix of a graph are related to the properties of the graph. This is the topic of spec- 
tral graph theory. A hypergraph can be described by a (0, 1)- symmetric tensor, which 
is called its adjacency tensor. Are the properties of the eigenvalues of the adjacency 
tensor of a hypergraph related to the properties of the hypergraph? Recently, several 
papers appeared in spectral hypergraph theory [91 [10l |2U SQl [51]. It is interesting that 
in [9j QUI El], H-eigenvalues are used, while in [lOj El], Z-eigenvalues are used. This 
shows that there are more to explore in this topic. We expect a comprehensive spectral 
hypergraph theory will emerge eventually. 

(4) The Relation with Nonnegative Tensor Factorization. Nonnegative tensor 
factorization is an important topic in the tensor decomposition community [22} El]. What 
is the relation between the spectral theory of nonnegative tensors with the nonnegative 
tensor factorization topic? Currently, it is still blank in the literature. 

C. Eigenvalues of General Tensors: Theory 

This group include six research topics. 

(1) The Number of E-Eigenvalues, the Number of Real Eigenvalues, E- 
Characteristic Polynomial. In [65J, Ni, Qi, Wang and Wang showed that the degree 
of the E-characteristic polynomial of an even order tensor is bounded by 

, (m-l) n -l 
d{m, n) = 



m — 2 

where m is the order and n is the dimension, m is even and m > 4. The definition 
of E-eigenvalues in [70J is not a strict generalization of eigenvalues of a square matrix, 
as it excludes the case that x T x = 0. In [llj, Cartwright and Sturmfels gave a precise 
definition by using equivalence classes, and showed the number of the equivalent classes 
is 

(m - 1)" - 1 

d[m,n) = , 

m — 2 

for any m > 2, echoing the result of [65] . 

In [H] , Chang, Pearson and Zhang unified the definitions of H-eigenvalues, Z-eigenvalues 
and D-eigenvalues, and showed that for a real n- dimensional mth order symmetric ten- 
sor, there are at least n H-/Z-/D-eigenvalues. Existence of real eigenvalues of real tensors 
was further studied in |105j . 

In [50], Li, Qi and Zhang studied the properties of E-characteristic polynomials. 

(2) Symmetric Hyperdeterminant or E-determinant. In [2l], Cooper and 
Dutle gave explicit expression of symmetric hyper-determinants of some special classes 
of hypergraphs for arbitrary n. This shows that for some special cases, symmetric hyper- 
determinants are computable. In [50], Li, Qi and Zhang showed that symmetric hyper- 
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determinants are invariant under orthogonal transformation. In |36j, symmetric hyper- 
determinants are called E-determinants. By the definition, det(*4.) = if and only if 
Ax m ~ l = has a nonzero solution. Hu, Huang, Ling and Qi [36] showed that if det(vA) 7^ 
0, then any nonhomogeneous polynomial system with A as its leading coefficient tensor 
has a solution. This shows that like determinants are closely related to the solution 
theory of linear systems, hyper-determinants are closely related to the solution theory of 
polynomial systems. More research is needed in this topic. 

(3) Singular Values, Symmetric Embedding and Eigenvalue Inclusion. Sin- 
gular values of rectangular tensors were studied in [18j [191 E0 EH1 HOlj 1111] . Symmetric 
embedding, which links singular values of rectangular values, and eigenvalues of square 
tensors, was studied in [20j |82] . 

In [5T1 152] . the authors give a number of eigenvalue inclusion theorems, including 
a Taussky-type boundary result and a Brauer eigenvalue inclusion theorem for tensors, 
and their applications. 

(4) Waring Decomposition Related Problems. A Waring decomposition of a 
(homogeneous) polynomial / is a minimal sum of powers of linear forms expressing /. 
Under certain conditions, such a decomposition is unique. In |67J, Oeding and Ottaviani 
discussed some algorithms to compute the Waring decomposition, which are linked to the 
equations of certain secant varieties and to eigenvectors of tensors. In particular they 
explicitly decomposed a cubic polynomial in three variables as the sum of five cubes 
(Sylvester Pentahedral Theorem). 

(5) Third-Order Tensors as Linear Operators. In [H[30l|45], special operations 
were introduced to regard third-order tensors as linear operations. With more compli- 
cated manipulations, this approach may also be extended to higher-order tensors. The 
cost of this approach is that the original tensor operations are given up. 

(6) Successive Symmetric Best Rank-One approximation, the Best Rank- 
One Approximation Ratio, and the Symmetric Rank Conjecture. Denote the 
space of all real mth-order n dimensional symmetric space as S m ^ n . Suppose that A G 
S m>n . Let ^4^°^ = A. In general, suppose that A^ G S m>n . Let be the Z-eigenvalue of 
A^ with the largest absolute value, and be the corresponding Z-eigenvector. Then 
we define 

A^=A^-\ k 

Wang and Qi [92] proved that 

-4=5>(^r, 0) 

and called this the successive symmetric best rank-one approximation of A. 

Qi [73] studied further the convergence of ([9]). Denote the largest absolute value of 
the Z-eigenvalues of A G S mtn as pz(A). When m = 2, pz(-A) is actually the spectral 
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radius, i.e., the 2-norm of A. In S m ,m Pz(A) is still a norm of A. Denote the Frobenius 
norm of A by Then 

App{S m , n ) : = min : A G S m A > 0. 

Qi [73] called it the best rank-one approximation ratio of 5" m n . Actually, we have 

< App(S m , n ) < 1. 

Suppose that A^ G S m>n . Let be the Z-eigenvalue of A^ with the largest absolute 
value, and be the corresponding Z-eigenvector. Then we have 

||^ (fe+1) ||^<||^ (fe) ||^[l-^p(5 m , n ) 2 ]. 

This not only proves (El), but also gives its convergence rate. 
Let 

1 

/4 



— m,n ^f, 



n 



m—l 



If m = 2k is even, then let A^" 1 '^ G and [i m ,n De defined by 

A M x m = (x T x) k (10) 



and 



fJ>m,n 



I M(m,n) II ' 

If m = 2/c + 1 is odd, then let ^4( m ' n ) e 5 m)n and /2 OT)n be defined by 



A^x m = (x T x) \J2 Xi ) ( n ) 

and 



,i=i 



{J"m,r 



||^(m,n) || ' 

The following upper and lower bounds for App(S m ^ n ) were established in 

Theorem 4. The value /i a positive lower bound for App(Ss jTl ). 

On the other hand, the value p, m ,n i> s an upper bound for App(S mtn ) for m = 2, 3, • • • . 
We have ^ 

— = K „ = App(S 2 , n ) = fi2,n = —j=, (12) 

'n — 2 ' n \/n 



and 



1 / 6 

- = Mo _ < App(S m , n ) < fx 3tn = \ — — (13) 
n — d ' n \ n + o 



A P p(S 4 , n ) < /2 4 , n = \l n 2 + 2n ( 14 ) 



9 



To improve this theorem, Qi [73] made a conjecture that pz(A) is equal to 

a(A) :=max{|Ac (1) ---a; (m) | : (x (i) ) T x (l) = 1, Vi = 1, • • • , m} . 

This conjecture was proved by Zhang, Ling and Qi |108j . This improved Theorem H] that 
fi is a positive lower bound for App(S m ,n) for all m. 

One referee of [108J pointed out the above result can be regarded the first step to 
prove the symmetric rank conjecture of Comon, Golub, Lim and Mourrain [23J. 

Let A be an mth-order n-dimensional tensor. We may have a rank-one tensor de- 
composition of A: 

A=J2a lX ^---x^ m \ 

i=l 

where x^'^ are n-dimensional vectors, ctj are numbers. The minimum value of r is called 
the rank of A. 

If A is symmetric, than we may have a symmetric rank-one tensor decomposition of 

A: 

A = jz^ (* ( T> 

where are n-dimensional vectors, «j are numbers. The minimum value of r is called 
the symmetric rank of A. 

It is conjectured by Comon, Golub, Lim and Mourrain [23] that for a symmetric 
tensor, its symmetric rank is the same as its rank. We call this conjecture the symmetric 
rank conjecture. 

D. Eigenvalues of General Tensors: Links with Optimization, Numerical 
Analysis and Geometry 

This group include four research topics. 

(1) Bi-quadratical Optimization and Spherical Optimization. Bi-quadratical 
optimization and spherical optimization were studied in [7J [33j [351 [591 EDI EH |85l |93| 
1107} \108\ 1110] . Bi-quadratical optimization is linked with the M-eigenvalue problem 
[?, ?,|22J[H3], while spherical optimization is linked with the Z-eigenvalue problem. They 
are NP-hard in general. The approach in [7J [33j |59j |60j |63j [851 H07[ 1108] is to use semi- 
definite programming to find a lower bound for the minimization problem. Symmetric 
spherical optimization was studied in |108[ 1110] . As stated in the last section, Zhang, 
Ling and Qi |108] proved that for A E S m>n , 

p z (A) = max {|Ar m | : x T x = l} 

= a(A) =max^\Ax {1) ■■■x (m) \ : (x^Y x® = 1, Vi = 1, • • • ,m}. 
This reveals some deep insight property of symmetric spherical optimization. 
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(2) Space Tensor Conic Programming. In the study of higher order diffusion 
tensor imaging, positive semi-definite tensor models arise [37J EOj 151] . where the tensor 
is positive definite and of dimension three. We call three dimensional tensors space 
tensors. Space tensors appear in physics and mechanics. They are real physical en- 
tities. The positive semi-definite space tensors of the same order form a convex cone. 
The optimization problem involving positive semi-definite space tensors is called the 
space tensor conic programming problem [79]. In [31], based on the analysis of the 
semismoothness properties of the maximum Z-eigenvalue function, a generalized Newton 
method is proposed to solve the space tensor conic linear programming problem. In |39j, 
a sequential SDP method is proposed to solve the space tensor conic linear programming 
problem. 

The problem to identify a space tensor is positive (semi-) definite or not is equivalent 
to finding the minimum Z-eigenvalue of an ra-th order tensor with dimension three. As we 
discussed in the basic theory section, this problem can be solved by elimination methods. 
It is not NP-hard. This approach requests to solving a one variable polynomial equation 
of degree generically iV = |(m + l)(m + 2). When m is big, the solution procedure is 
unstable. Hence, new algorithms are needed and this research topic is interesting. 

(3) Characterizing the Limiting Behaviour of Newton's Method. Dupont 
and Scott [26] studied the angular orientations of convergent iterates generated by New- 
ton's method in multiple space dimensions. They showed that the Newton iteration can 
be interpreted as a fixed-point algorithm for solving a tenor eigenproblem. They gave 
an extensive computational analysis of this tensor eigenproblem in two dimensions. In 
a large fraction of cases, the tensor eigenproblem has a discrete number of solutions to 
which the Newton directions converge quickly, but there is also a large fraction of cases 
in which the behavior is more complicated. Dupont and Scott contrasted the angular 
orientations of iterates generated by Newton's method with the corresponding directions 
of the continuous Newton algorithm. 

(4) Geometry related Problems. Qi [71] studied the orthogonal classification 
problem of real hypersurfaces given by an equation of the form 

S={xeM n :f(x)=Ax m = c}, 

where A is a symmetric tensor, through the rank, the Z-eigenvalues and the asymptotic 
directions. 

Balan [TJ |2J El E] showed that there is a close relationship between Finsler metrics 
of Berwald-Moor type, and the spectral theory of tensors. Balan used the results of 
[70| [71] to analyze spectral properties of m-root models, geometric relevance of spectral 
equations, degeneracy sets, asymptotic rays, base indices and based rank one approxi- 
mation of various Berwald-Moor tensors, Chernov tensors and Bogoslovski tensors. A 
further investigation on the geometrical meanings of eigenvalues of higher order tensors 
is needed. 
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E. Eigenvalues of General Tensors: Links with Image Science, Solid Me- 
chanics and Quantum Physics 

This group include four research topics. 

(1) Higher Order Diffusion Tensor Imaging. Diffusion tensor imaging (DTI) is 
the most popular magnetic resonance imaging model [5]. In the DTI model, eigenvalues 
of the diffusion tensors play a fundamental role as they are the main invariants of these 
tensors, and thus can be used in the measures of DTI [5]. However, DTI is known to 
have a limited capability in resolving multiple fibre orientations within one voxel. This 
is mainly because the probability density function for random spin displacement is non- 
Gaussian in the confining environment of biological tissues and, thus, the modeling of 
self-diffusion by a second order tensor breaks down. Hence, researchers presented various 
Higher Order Diffusion Tensor Imaging models to overcome this problem [HJ [621 
EEJ EU ES] • For the measures of these higher order diffusion tensor imaging models, the 
spectral theory of tensors naturally play an important roles [6j EH [321 [371 H21 E31 UM 
[781 1801 EH 1109] . In particular, Qi, Wang and Wu [78] proposed D-eigenvalues to study 
the diffusion kurtosis model [HI [62]; Bloy and Verma [6] studied the underlying fiber 
directions from the diffusion orientation distribution function by using Z-eigenvalues; Qi, 
Yu and Wu [81] . Hu, Huang, Ni and Qi [37J, and Qi, Yu and Xu [80] studied positive 
semi-definite higher order diffusion tensor imaging models. 

(2) Image Authenticity Verification Problem. In [100] . Zhang, Zhou and Peng 
generalized the definition of D-eigenvalues, and introduced the gradient skewness tensor 
which involves a three-order tensor derived from the skewness statistic of gradient images. 
They found out that the skewness value of oriented gradients of an image can measure the 
directional characteristic of illumination, and the local illumination detection problem 
for an image can be abstracted as solving the largest D-eigenvalue of gradient skew- 
ness tensors. Their method presented excellent results in a class of image authenticity 
verification problem, which is to distinguish real and flat objects in a photograph. 

(3) Elasticity Tensors in Solid Mechanics. In solid mechanics, the elasticity 
tensor A = (a^n) is partially symmetric in the sense that for any i,j,k,l, we have 
dijki = O'kja = a ukj- We say that they are strongly elliptic if and only if 

n 

f(x, y) = Axyxy = ^ <kjhffliyphVi > 0, 

i,j,k,l=l 

for all unit vectors x, y G 3ft", n = 2 or 3. For an isotropic material, some inequalities 
have been established to judge the strong ellipticity [3E1 [91]. In [311 Dai, Han and 
Qi studied conditions for strong ellipticity and introduced M-eigenvalues for the ellip- 
ticity tensor A. They showed that M-eigenvalues always exist and the strong ellipticity 
condition holds if and only if the smallest M-eigenvalue of the elasticity tensor is positive. 

(4) Quantum Entanglement Problem. The entanglement problem is to de- 
termine whether a quantum state is separable or inseparable (entangled), or to check 
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whether an run x ran symmetric matrix A >z can be decomposed as a convex combina- 
tion of tensor products of n and m dimensional vectors. It has fundamental importance 
in quantum science and has attracted much attention since the pioneer work of Einstein, 
Podolsky and Rosen [27] and Schrodinger [83]. The typical optimization problem in 
quantum entanglement problem is like the following [34J: 

max { • ■ ■ x (m) | : \x® \ — 1, for i — 1, • ■ • , m} , 

where A is an mth-order (ni, • ■ ■ , ra m )-dimensional complex tensor, x^' E C rH for i = 
1, ■ • • , m. The singular value form of this optimization problem is 

Ax w ■ ■■x {i ~ 1) x {i+1) ■■ -x {m) = Xx {i) , \x®\ = 1, for i = 1,-- • ,m. 

The optimization problem and the singular value problem involve complex and conjugate 
numbers. They are different but have some similarities with what we studied before. 
With their strong physics background, this topic will be a new fertile land of our research. 

The International Conference on The Spectral Theory of Tensors will be held at 
Chern Research Institute, Nankai University, Tianjin, China, during May 30 - June 2, 
2012. Its website is: 

http://www.nim.nankai.edu.cn/activites/conferences/hy20120530/index.htmj 
Two special issues on the Spectral Theory of Tensors are being edited [5HI I5B] . 
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